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Abstract 

o 

■^j- ■ The competition between commensurate and incommensurate spin-density- 

wave phases in the infinite-dimensional single-band Hubbard model is exam- 
ined with quantum Monte Carlo simulation and strong and weak coupling 
approximations. Quantum fluctuations modify the weak-coupling phase dia- 

o" 

O ■ gram by factors of order unity and produce remarkable agreement with the 



quantum Monte Carlo data, but strong-coupling theories (that map onto effec- 



C$ • tive Falicov-Kimball models) display pathological behavior. The single-band 

model can be used to describe much of the experimental data in Cr and its 
dilute alloys with V and Mn. 
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Introduction. Spin-density-wave (SDW) order, in which the modulation wavevector of 
the SDW is incommensurate with the underlying lattice, is one of the most fascinating 
ordered states found in nature. Incommensurate magnetism occurs in both metallic and 
insulating phases and on both frustrated and unfrustrated lattices. In general, incommen- 
surate magnetic order may be driven by either frustration or by Fermi surface nesting with 
a wavevector that lies away from commensurate wavevectors. It is important to understand 
which process plays a more important role, and to understand how many-body effects mod- 
ify the stability of incommensurate phases. Here the effect of nesting is examined on an 
unfrustrated lattice with strongly correlated electrons. The resulting phase diagram is then 
compared to approximate results in the weak- and strong-coupling limits (Figs. 2 and 3). 
Finally, our theoretical results are compared to those found in Cr and its dilute alloys. 

Elemental Cr is a paradigm for an antiferromagnetic (AF) metal J1J with incommensurate 
order driven by Fermi-surface nesting. The lattice structure of Cr is an unfrustrated body- 
centered-cubic structure which may be modeled by a Hubbard model || near half filling 
with moderate electron-electron correlations. Adding electrons to Cr (by alloying with Mn) 
rapidly makes the magnetic order commensurate with the lattice, whereas removing electrons 
from the system (by alloying with V) rapidly increases the incommensuration and decreases 
the magnetic transition temperature, eventually to zero ||. 

Heretofore, incommensurate magnetic order has mainly been examined within the 
Hartree-Fock (HF) (weak-coupling) approximation that neglects quantum fluctuations. 
Penn |4j] found incommensurate order in the three-dimensional Hubbard model, and 
Schulz [[J found evidence for incommensurate phases on a square lattice. 

Formalism. In this contribution the magnetic phase diagram of the single-band Hub- 
bard model is investigated in the limit of infinite dimensions ||. This limit is useful, because 
it has been shown to contain most of the physics expected of three-dimensional Hubbard 
models, and the many-body effects can be treated essentially exactly with the quantum 
Monte Carlo (QMC) techniques of Hirsch and Fye j?|. This allows us to demonstrate the 
existence of incommensurate order at finite temperatures in a model that only includes 



Fermi-surface nesting effects and electron-electron correlations. We find incommensurate 
magnetism to be a ubiquitous property of this model. 

The Hubbard model || is described by the following Hamiltonian: 
t* 11 
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where c\ (q) is a creation (destruction) operator for an electron at site i with spin a. The 
hopping matrix elements connect nearest neighbors on a hypercubic lattice in ^-dimensions; 
its magnitude is written as t = t* /2yd [to have a well-defined limit in infinite dimensions 
(d — > oo)]. All energies are expressed in units of the rescaled hopping matrix element t*. 
The Coulomb repulsion is represented by U and the chemical potential by ji. 

In the limit of infinite dimensions the local approximation becomes exact. The electronic 
Green's function G(iu n ) = G n is represented by an integral over the noninteracting density 
of states p(y) = exp(—y 2 )/^/Tr 

G n = dy~ = Foo^Un + fj, - E„) , (2) 

J-oo iu n - fi - L n - y 
with S n = S(zu; n ) the electronic self energy. The magnetic susceptibility satisfies Dyson's 
equation 

Xmn(q) = Xm(q)^mn ~ T ^ Xmi^ mrXrn(q) , (3) 

r 

with T mr = T(iu m , iu r ) the (local) irreducible vertex function in the SDW channel. The 
bare susceptibility x^(q) is defined by 



X ^(q) = -T5:G m (k)G m (k + q) 
k 
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The bare susceptibility only depends on the scalar parameter X(q) = J3f =1 cos qj/d, which 
defines an equivalence class of wavevectors in the infinite-dimensional Brillouin zone ||. 
X(q) can be parametrized by the line that extends along the Brillouin zone diagonal from 
the zone center (X = 1) to the zone corner (X = —1). The self energy and irreducible vertex 



function are extracted from the self-consistent QMC simulations as described previously || 
(for each value of U, we employ a constant small value of At for the spacing of the imaginary 
time slices in the Trotter breakup). 

It is important to qualify the region of validity of different approximation techniques 
by comparing them with the essentially exact QMC results. Here we compare both weak- 
coupling and strong-coupling perturbation theories. 

In the weak coupling limit a renormalized Hartree-Fock approach |10] is employed. The 
Neel temperature is determined by the Stoner criterion 

i = X°(X,T)=T$>°(X) , (5) 

U n 

where the bare particle-hole susceptibility XnPO * s calculated with noninteracting Green's 
functions G° [E n = in Eq. (|)]. This HF transition temperature is reduced by factors 



of order three due to quantum fluctuations |pL0|| , even in the limit U/t* — > 0. Quantum 



fluctuations modify the Stoner criterion to WT 



1 = X (X,T)- X (X = 0,T) = X °(X,T) + T£(G°) 2 , (6) 

U n 

in the limit U/t* — » 0. These fluctuations initially reduce T c by the factor 



exp[— x°(0, T c )/p(n)] in the weak-coupling limit ||T0| . |TT . 

In the strong-coupling limit the Hubbard model can be mapped onto a Falicov-Kimball 
model fl2,13|. This mapping is exact in the limit U/t* — > oo (which reproduces the atomic 



limit) for the self energy but is not exact for the irreducible vertex functions in this limit. As 
a result, the strong-coupling theories display pathological behaviors. More explicitly, these 
approximations assume that the down-spin particles form a static background when the up- 
spin particles move and vice versa; this system is then described by Falicov-Kimball models 
fT4| for both the spin-up and spin-down electrons that are self-consistently coupled together. 
Two different self-consistent coupling schemes have been proposed so far |T2| . |T3|| . Janis 
and Vollhardt's approximation JT2J underestimates the SDW susceptibility at half filling 



which strongly suppresses T c and does not reproduce the Heisenberg limit of T c ps t* 2 /2U. 



Li and d'Ambrumenil's approximation Jl3| is correct for large U at half filling, but has 
the pathological behavior of predicting ferromagnetism away from half filling (with a finite 
transition temperature in the U/t* — * oo limit) because of segregation in the effective Falicov- 
Kimball model. This latter pathology occurs because the zero temperature occupation 
number of the static particles is 0, 0.5, or 1, and segregation occurs in the Falicov-Kimball 
model whenever the static particle concentration is 0.5, and does not equal the mobile 



particle concentration |15| . We deal with this pathology, by only considering ordered states 
with X < 0. 

Results. To determine the magnetic transition temperatures, we calculated the mag- 
netic susceptibility for all X(q) in the Brillouin zone. As shown in Fig. 1, the susceptibility 
always displayed a maximum at a distinct value of X. The transition temperature was then 
inferred from interpolation (or extrapolation) of the peak inverse susceptibility, as shown 
in the inset to Fig. 1. At half filling || the Hubbard model in infinite dimensions has a 
transition to a commensurate AF state (X = —1) at the Neel temperature T^. T^ increases 
monotonically from zero at U = to 0.14t* at U « 3t*. When U is increased further, 
T/v decreases monotonically to zero at U = oo. As the system is doped away from half 
filling, the Neel temperature drops until a critical filling is reached where the commensurate 
SDW becomes incommensurate. This is shown in Figs. 2a and 2b for the weak-coupling and 
strong-coupling results respectively. As the system is doped further away from half filling, 
the wavevector of the ordered phase changes continuously with the electron concentration 
until T c drops to zero at the incommensurate-paramagnetic phase boundary. The shape of 
the magnetic phase boundary changes continuously from a BCS-like curve (as a function 
of doping) at weak coupling to an almost linear curve at strong coupling, with the same 
crossover region (U « 3t*) as found for T/v at half filling. 

Phase diagrams in the weak-coupling regime (U/t* < 3) have been obtained from both 
QMC simulations and from the theory of Eq. (^). They are plotted in Fig. 2(a) for four 
values of U/t* (U/t* = 1, 1.5, 2, 3). Let X max denote the largest value of the scalar parameter 
X with which incommensurate order is found for each value of U. In the QMC simulations, 



we find that both X max increases, and that the ratio of the transition temperature at the 
commensurate- incommensurate phase boundary Tj to the Neel temperature at half filling T^ 
decreases, as the coupling strength increases. However, when the modified Stoner criterion 
is used, we find that the transition temperature curves scale with the coupling strength and 
maintain the same approximate shape. This implies that X max will increase, while the ratio 
Ti/Tn remains constant, Tj/Tjv ~ 0.57. 

QMC results in the strong-coupling regime (U/t* > 3), along with Li and d'Ambrumenil's 
approximation for T c are plotted in Fig. 2(b) for U/t* = 3,4,5,7. The approximate results 
are generated with the restriction that only ordered states with X < are considered 
(which suppresses the ferromagnetism due to phase separation). Under this assumption, the 
approximate strong-coupling theory predicts no incommensurate order (near X = — 1), and 
the transition temperature curves also maintain the same shape as the coupling strength 
changes. Accurate simulations at very large values of U are not possible with the QMC. 
Thus, we are unable to determine whether X max continues to increase in the strong-coupling 
regime, nor are we able to determine what happens to Tj/Tjy. Of course, the strong-coupling 
approximation predicts both X max and Tj/T^ are constants, since it does not display any 
incommensurate order. 

Comparison of both weak-coupling and strong-coupling approximations with the QMC 
solutions shows that the approximate methods are unable to reproduce the qualitative 
change in shape of the finite-temperature phase diagrams as a function of U. 

The "phase diagram" which indicates the commensurate-incommensurate phase bound- 
ary (occurring at T = Tj) and the incommensurate-paramagnetic phase boundary (occurring 
at T = 0) is presented in Fig. 3. The thin (thick) solid lines denote the commensurate- 
incommensurate boundary for the Stoner (modified Stoner) criterion; the thin (thick) 
dashed lines plot the corresponding paramagnetic phase boundary. The dotted line is the 
commensurate-paramagnetic phase boundary in the strong-coupling theory, and the dots are 
the QMC results. The quantum fluctuations strongly renormalize the HF phase boundary 
to produce good agreement with the QMC. The value of X max = occurs at U — oo when 
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the modified Stoner criterion is used, whereas X max increases to the ferromagnetic point 
(X = 1) in the Stoner theory (ferromagnetism does occur at large values of U and small 
electron concentration when the modified Stoner criterion is used). Note that the changes 
to the phase diagram brought about by the fluctuations are simple: the critical value of U 
is shifted by 1/U = 1/Ujjf ~ X°(0> T c ), which pushes some of the anomalous behavior of the 
HF approximation beyond U = oo. The QMC calculations show that many-body effects 
actually favor the formation of incommensurate phases, since the region of stability of in- 
commensurate magnetism is larger than what the weak-coupling or strong-coupling theories 
would predict. 

Comparison to Experiment. Finally, it is of interest to compare our results with what 
is known about magnetism in elemental Cr. Electronic band structure calculations [1] show 
that the <i-electron concentration for Cr is 4.6/atom which is close to a half-filled band. 
Doping with Mn adds an electron to the <i-bands, and doping with V removes an electron. 
The commensurate-incommensurate phase boundary lies at a doping of 0.3% Mn and the 
paramagnetic phase boundary at a doping of 3.5% V |j]||. Since the density of states for Cr 
is peaked near the band edges, rather than the band center, it is difficult to map directly onto 
the Gaussian density of states (of the single-band model in infinite dimensions). Instead, 
we compare the ratio of the incommensurate ordering wavevector to the commensurate 
wavevector in order to estimate the magnitude of the Coulomb interaction for the effective 
single-band model. The smallest value for the ratio of the incommensurate wavevector to 
the commensurate wavevector is 0.92 for Cr |]|| implying X max = cos(0.927r) = —0.97 for 
the single-band model. The approximate value of U is then estimated to be U ss 1.9t* 
(U ~ 1.5t*) for the weak-coupling theory (QMC). (We also found that the incommensurate 
wavevector changes very rapidly with doping near the commensurate-incommensurate phase 
boundary which is reminiscent of the first-order jump in the wavevector that is seen in Cr.) 
The ratio of Tj (325i^) to T/v (700.fQ is 0.46 which is smaller than the weak-coupling value 
of 0.57 and is consistent with the assignment of a small value to U/t*. The Neel temperature 
at half filling is approximately 0.096 t* (0.086 t*) in the weak-coupling theory (QMC) which 



yields an effective bandwidth W = At* = 2.4 eV (2.7 eV) for the single-band model. This 
is a reasonable number since Cr has a bandwidth of 6.8 eV and the pileup of the density of 
states at the band edges JI[, implies that the bandwidth for the effective single-band model 
must be larger than the naive approximation of one-fifth of the total bandwidth. 

Conclusion. We have shown that incommensurate SDW order exists in the infinite- 
dimensional Hubbard model as one dopes away from half filling. A simple modification 
(due to quantum fluctuations) of the usual Stoner criterion produces good agreement with 
the weak-coupling QMC results and is easy to implement in arbitrary dimensions (the q- 
dependent susceptibility is reduced by the local susceptibility before applying the Stoner 
criterion). Since these quantum fluctuations produce large renormalizations of the magnetic 
phase boundaries, it is worthwhile to repeat previous Hartree-Fock calculations ||||, and 
employ the modified Stoner criterion, to more accuratly determine the phase diagram in the 
weak-coupling regime. We have also found that strong-coupling theories (which map onto 
effective Falicov-Kimball models) display pathological behavior in the magnetic transition 
temperature, because the (approximate) irreducible vertex functions do not reproduce the 
atomic limit when U/t* — > oo. Finally, we have shown that much of the behavior found in 
Cr and dilute Cr alloys can be described by an effective single-band Hubbard model that 
does not include any of the microscopic details of the bandstructure. 
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FIGURES 
FIG. 1. The magnetic susceptibility for all X(q) at various temperatures when U = 4 and 

p e = 0.825. The susceptibility displays a peak at X ~ —0.90. As shown in the inset, the transition 

temperature (T c = 0.0148i*) was inferred from extrapolation of the peak inverse susceptibility. 
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FIG. 2. Phase diagram of the Hubbard model in the (a) weak-coupling regime 
(U/t* = 1, 1.5, 2, 3) and (b) strong-coupling regime (U/t* = 3, 4, 5, 7). The solid (open) dots denote 
the transition temperature to a commensurate (incommensurate) SDW phase as determined by a 
QMC calculation. The solid (dotted) lines denote the transition temperature to a commensurate 
(incommensurate) SDW phase using the modified Stoner criterion in Figure 2 (a) and using Li and 
d'Ambrumenil's approximation in Figure 2 (b). The dashed lines are a fit of the QMC data to the 
form T c /t* = a(x — b) c . The exponent c increases with increasing U. 
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FIG. 3. "Phase diagram" for the Hubbard model as a function of electron concentration and 
U. The thin (thick) solid lines denote the commensurate-incommensurate phase boundary for the 
Stoner criterion (modified Stoner criterion); the thin (thick) dashed lines are the corresponding re- 
sults for the incommensurate-paramagnetic phase boundary. The dotted line is the strong-coupling 
approximation for the commensurate-paramagnetic phase boundary. The solid (open) dots denote 
the QMC solutions that display commensurate (incommensurate) SDW order. 
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